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Abstract 

The World Trade Web (WTW), the network defined by the international im- 
port/export trade relationships, has been recently shown to display some impor- 
tant topological properties which are tightly related to the Gross Domestic Product 
of world countries. While our previous analysis focused on the static, undirected 
version of the WTW, here we address its full evolving, directed description. This 
is accomplished by exploiting the peculiar reciprocity structure of the WTW to 
recover the directed nature of international trade channels, and by studying the 
temporal dependence of the parameters describing the WTW topology. 
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1 Introduction 



The internal structure of many biological, social and economic systems dis- 
plays complex topological properties such as clustering, a scale-free degree 
distribution and the presence of degree correlations, which are not reproduced 
by simple random graph models [1]. The global trade activity is one example of 
a large-scale system whose internal structure can be represented and modeled 
by a graph [2,3]. In the most general case, from publicly available annual trade 
data [4] it is possible to define the evolving World Trade Web (WTW in the 
following) , where each world country is represented by a vertex and the flow of 
money between two trading countries is represented by a directed link between 
them: if % imports some good from country j during the year t, then a link from 
i to j is drawn in the t-th snapshot of the graph, corresponding to a nonzero 
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entry a,ji(t) = 1 in the corresponding adjacency matrix. Otherwise, no link is 
drawn from % to j and %i(£) = 0. Note that the direction of the link always 
follows that of the wealth flow, since exported (imported) goods correspond 
to wealth flowing in (out). In such a description, if N(t) denotes the number 
of world countries during year t, the in-degree k™(t) = J2f=i a ij(t) and the 
out-degree k° ut (t) = J2j=i a ji{t) of a country % correspond to the number of 
countries importing from and exporting to i respectively. 

The structure of the WTW continuously changes in time due to variations 
in the number N(t) of world countries and to the continuous rearrangement 
of its connections, and hence of (^ (i). While previous studies focused on the 
undirected description of a single snapshot of the WTW [2,3], in the present 
work we study the properties of the WTW as a directed and evolving network 
and we confirm that its topology always displays a peculiar dependence on 
the Gross Domestic Product of world countries, in excellent accordance with 
the so-called fitness model [5]. We first review the previous results and then 
extend them by reporting the temporal evolution of the relevant parameters. 
Our analysis is based on a comprehensive dataset [4] reporting the annual 
trade activity between all world countries together with the annual values of 
their total Gross Domestic Product during the period 1950-1996. 



2 Topology of the WTW 



Two recent papers document some interesting static topological properties 
of the WTW viewed as an undirected network [2,3]. In such a representa- 
tion, each pair of vertices i,j is considered connected if % and j trade in at 
least one direction. Before reviewing these results, we first give the explicit 
expressions linking the directed and undirected versions of a graph. From a 
graph-theoretic point of view, if is the adjacency matrix of a directed 
graph, then the undirected version of the same graph is described by the sym- 
metric adjacency matrix 6^ = + — a^a^. Correspondingly, the quanti- 
ties k™ = J2j a ij an d k° ut = J2j a ji computed on the directed graph and the 
'undirected' degree k{ = J2j % computed on its undirected version are related 
through the following equation: 

h = kT + kT l - k~ (i) 



where we have introduced the reciprocal degree [6] kp = J2j o-ijO-ji, defined as 
the number of neighbours of % with double (reciprocal) connections to and 
from %. While in principle there is an obvious loss of information in using the 
undirected version of a directed graph, in a recent paper [6] we showed that for 
the particular case of the WTW there are two important empirical properties 
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Fig. 1. Properties of the WTW. a) Reciprocal (k*~*) versus total (k T ) degree for 
various timesteps. b) Time evolution of the reciprocity r(t). c) Time evolution of 
the number of world countries N(t). 

that allow to recover the directed character from its undirected description. 
Firstly, one always observes that approximately kl n (t) w k° ut (t) Vi, t. Sec- 
ondly (see fig. la), a linear relation holds between k£*(t) and the total degree 
kf(t)=kt(t) + kr(t) [6]: 



(2) 



where r(t) is the reciprocity, defined as the ratio of the number L*~*(t) of links 
pointing in both directions to the total number of links L(t): 



r(t) = 



L(t) 



(3) 



The annual values of r(t) for the 1950-1996 time interval are shown in fig. lb. 
The above results allow to obtain the properties of the WTW, viewed as a 
directed graph, from its undirected version. For instance, eqs. (1,2) imply 



kTit) « kr\t) 



kj{t) „ h(t) 



2 - r(t) 



(4) 



and a relation between L and the number L u of links in the undirected network: 

1 N(t) 9 

E Ht) 



N(t) 

L(t) = y ki n (t) = 



2 - r(t) 



L u (t) 



(5) 



As a consequence, in the following we can simply address the topology of 
the undirected WTW and then recover, at any timestep t, its full directed 
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description through the value of r(t). For simplicity, we now review the results 
corresponding to a single snapshot (the year 1995) of the undirected WTW 
[3]. An important result that we report here is that the same qualitative trends 
are observed for each year in the database (1950-1996) and can be described 
quantitatively by plotting the temporal evolution of the relevant parameters. 
For instance, in fig.lc we show the evolution of the number of countries N(t) 
during such time interval. For the 1995 data, the number of countries equals 
N = 191. 

One first-order topological property is the degree distribution P{k). In ref. [3] 
we showed that, contrary to previous results [2] suggesting a power- law form 
of the distribution with exponent —1.6, the scale-free behaviour is suppressed 
by a sharp cut-off, which is actually restricted to a narrow region (see fig. 2a, 
where a power law with exponent —1.6 is shown as a comparison). 



A second-order topological property, related to the form of the two-vertices 
joint degree distribution P(k,k'), is obtained by plotting the average nearest 
neighbour degree (ANND), defined as K™ n = -^-J2j^ijkj, versus the degree ki 
(see fig. 2b). The decreasing trend shown in the figure clearly signals that the 
degrees of neighbouring vertices are anticorrelated, or in other words that the 
WTW is a disassortative network [7]. From an economic point of view, this 
means that countries with many trade partners are on average connected to 
countries with few partners [2,3]. 



Finally, the third-order correlation structure (related to the three- vertices joint 
degree distribution) can be investigated by plotting the clustering coefficient 
Ci = 

ki{ki-i) ^jfr ^k^i,j hjbjkbki versus the degree h. As shown in fig.2c, the 
trend is again decreasing, meaning that the partners of well connected coun- 
tries are less interconnected than those of poorly connected ones [2,3], a prop- 
erty sometimes referred to as hierarchy. 
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Fig. 2. Topological properties of the 1995 snapshot of the undirected WTW. a) Cu- 
mulative degree distribution P(k) and comparison with a power law with exponent 
—1.6 (dashed line), b) Plot of Kf n versus ki. c) Plot of Cj versus ki. 
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3 The fitness network model and the role of the GDP 



A very intriguing result [3] is that the above topological properties of the 
WTW are in excellent agreement with the predictions of the fitness network 
model [5], assuming that the probability that two vertices are connected is a 
function of some hidden variable or fitness associated to each of them (again, 
we are considering the WTW as an undirected network). In other words, 
if each vertex % is assigned a fitness value Xi, the model assumes that the 
probability pij that i and j are connected is a function p(xi,Xj) of Xi and 
Xj alone. All the expected topological properties therefore depend only on 
the functional form of p(x,y) and on the fitness distribution p(x). For the 
WTW, we showed that the Gross Domestic Product (GDP in the following) 
of a country can be successfully identified with the fitness of the corresponding 
vertex [3]. Extending the formalism to the evolving case, if Wi(t) denotes the 
GDP of country i at time t we can introduce the adimensional variable 

x .( t ) = = ^ (6) 



and interpret it as the fitness of vertex %. Then, once the form of p{x,y) is 
chosen, all the expected topological properties predicted by the model can be 
compared with the empirical ones. If at every timestep t the GDP is regarded 
as (proportional to) the potential number of trade connections that a country 
can develop, then the natural choice for p t (x,y) reads [8] 

pA^t),xM = 1 + mXi{t)Xj{t) (7) 



where S(t) > is the only free parameter of the model. The above form ensures 
that, at a given timestep, different realizations of the network with the same 
degree sequence {ki(t)}^Ji are equiprobable [8]. The parameter 5(t) can be 
fixed by requiring that the expected number of links in the undirected WTW 

1 N(t) 

^0 = «EEftfc(*W0] (8) 

»=1 3+i 



equals the observed one L u (t). After that, the topology is completely specified: 
the expected degree distribution is obtained by computing the expected degree 

Ht) = Y,Pt[xi(t),x j (t)] (9) 

3+i 
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while the expected average nearest neighbour degree is 



Hj^i Hk^j Pt [Xj(t) , Xj (t)]p t [Xj (t) , x k (t)] 



(10) 



and the expected clustering coefficient is 



Ci{t) 



Y.j+i Ek^j,i Pt [Xj (t) , Xj (t)]p t [Xj (t) , X k (t)]p t [x k (t) , Xj (t)] 

Ht)[h(t) - 1] 



(11) 



For instance, for the year 1995 we found that 5(1995) = 78.6 in order to have 
L u (1995) = L"(1995) [3]. The expected properties given by eqs.(9,10,ll) ob- 
tained with this choice are shown in fig.2 as solid lines superimposed to the 
empirical points. The accordance is excellent, indicating that the model indeed 
captures the basic aspects of the WTW topology. The comprehensive analysis 
of the whole data set that we are presenting here confirms that the same is 
true for each snapshot of the network. We can describe the time evolution 
in a compact way by studying the temporal change of the statistical GDP 
distribution and of the parameter S(t) (see figs.3a-b). More specifically, note 
that the discrete sums J2i on the r.h.s. of eqs. (9,10,11) can be replaced [5] by 
integrals of the form N(t) J dxp(x, t) where p(x, t) denotes the fitness distri- 
bution at time t. Therefore, besides S(t), the model predictions only depend 
on N(t) and p(x, t) and the evolution of these three quantities completely 
specifies the evolution of the network (together with r(t) if one is interested in 
the directed character of the WTW). In fig. 3a we plot the fitness distribution 
in its cumulative form p > (x,t) = p(y,t)dy. We find that the tail of the 
distribution always collapses to a time-independent Pareto tail with exponent 
-1, corresponding to a power-law p(x,t) pa p(x) oc x~ 2 . The evolution of S(t) 
during the 1950-1996 time interval is instead shown in fig. 3b. 



4 Discussion and conclusions 

A very important goal would be to relate the time variation of the parame- 
ters to the minimum possible number of external, unpredictable factors. For 
instance, the variation of N(t) is of course due to complicated geopolitical 
reasons and should therefore be taken as an input information. By contrast, 
the seemingly irregular behaviour of 5(t) can be probably explained in terms 
of more regular quantities. For instance, we recall that S(t) determines the 
number of links L u {t) = L(t) and is therefore related to the link density 
z (t) = N(t)[N(t)-i] ( rat i° °f observed to possible links). Interestingly, we find 
that z(t) displays relatively small fluctuations about the value 0.4 (see fig. 3c), 
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a result related to the way new links form in the network (for instance, when 
a country 'splits' in two or more countries, these can 'inherit' the links of 
the initial country). The form and evolution of p(x) has to be traced back to 
the dynamics of the GDP that takes place on the network. Models of wealth 
dynamics on complex networks [9,10,11,12] highlight how the topology deter- 
mines the form of the wealth distribution. The results discussed here show 
that the reverse is also true, and therefore lead us to the intriguing picture of 
a countinuous feedback between dynamics and topology, a scenario that inti- 
mately relates the results of network theory to those of economic modelling. 
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Fig. 3. a) Cumulative fitness distribution for four snapshots of the WTW (points) 
and comparison with a power-law with exponent -1 (solid line), b) Time dependence 
of the parameter 5(t) of the model, c) Time dependence of the link density z{t). 
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